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, We present an improved method of handling off-sheh effects in deep inelastic 

P\ _ 

nuclear scattering. With a firm understanding of the effects of the nuclear 

^ ', wavefunction, including these off-shell corrections as well as binding effects 

(D ! 

r~| , and nucleon-nucleon correlations, we can begin to examine the role of QCD 



in nuclei through an analysis of the moments of the nuclear inelastic structure 



^ \ function Our analysis is aimed on one side at extracting the depen- 

dence of the moments of the nucleon structure function by using the recent 
high X world Iron data and by properly subtracting nuclear eff^ects from the 
perturbative contribution. On the other side we compare quantitatively the 
behavior of the extracted moments with a simple 0{\/Q'^) phenomenological 
form and we determine the mass term for this parametrization. 
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Typeset using REVT^ 



The discovery of nuclear effects in Deep Inelastic Scattering (DIS) ultimately represents 
the fact that the quark and gluon structure of a free nucleon is different from that of a 
bound one. Nuclear effects have been found systematically in a number of high energy 
processes using both leptonic and hadronic probes (see e.g. and for recent reviews 
and references). Results are most often presented in the form of a ratio, Ra, of the nuclear 
cross section per nucleon to that of the free nucleon, where Ra = 1 within error bars means 
absence of nuclear effects. 

In this Letter we concentrate on the dependence of DIS in the intermediate and high 
Bjorken x regions (x > 0.3). Our aim is twofold: 1. On one side we investigate in detail the 
struck particle's off-shell effects which, within a class of models based on conventional degrees 
of freedom, can be considered as the major source of modifications of the dependence 
of the bound nucleon's structure function. In this context we make the observation that 
the proper framework to investigate off-shell effects in DIS from nuclei has to be sought 
within a relativistic theory and with a definite prescription for going beyond the Impulse 
Approximation (lA). As a side result, we show how one needs to be very careful about using 
ad hoc prescriptions. 2. On the other side we construct the quantities needed in a QCD 
analysis of the data, namely the (normalized) moments of the nuclear structure functions. 
We show that the off-shell effects investigated in the first part of the paper have a very 
small impact on the moments or that, in other words, nuclear effects can be factored out 
almost exactly. This allows us to extract the coefficients of the Higher Twist (HT) terms 
of the nucleon structure functions more accurately than previous determinations ( and 
references therein). 

1. While for the free nucleon it has been possible to address this question entirely within 
QCD a natural complication arises for the nucleus 0. On a phenomenological level 

we see that there are now two types of constituents, both the struck quark and the parent 
nuclear constituent, e.g. a nucleon, that can be off-shell. As we will show, the quark's tran- 
verse momentum and its virtuality are related differently in a free and in a bound nucleon, 
respectively and this affects their relative structure functions which, within a parton model 
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depend on the transverse degrees of freedom (or on the virtuahty) through the particles 
propagators. A careful treatment of off-shell effects along these lines, turns out to be crucial 
and in fact ad hoc procedures, such as the popular DeForest prescription |^ , fail consistency 
in the DIS region, e.g. by predicting a so far unobserved pattern of x-scaling violations. 

One can describe the DIS process for nuclear targets with a "double" lA, according to 
which the virtual photon scatters off a quark that in turn belongs to a nuclear constituent 
(e.g. a nucleon). lA is expected to work at the values of interest (Q^ > 2GeV^) and 
at X not in the proximity of and 1. By using a description in which the struck quark 
and nucleon transverse degrees of freedom are decoupled (see e.g. Ref. [0)' nuclear 
structure function can be written as a convolution formula in the longitudinal variables 

Fi{x)= jyzfA{z)F^{^^^, (1) 

where z = {kQ — k^)/Mv is the light-cone (LC) momentum fraction carried by the nucleon, 
being the bound nucleon's energy which is obtained from energy conservation at the nuclear 
vertex: k^ = Ma — \J Mj^^ + k"^ and ^3 = k ■ q. Q fA{z) is the nucleon's LC momentum 
fraction distribution, and F2 is the nucleon structure function. Calculations based on 
Eq.(|l]) reproduce the slope of Ra at intermediate x (0.3 < x < 0.55), provided a description 
of nuclear dynamics which properly accounts for Nucleon-Nucleon (NN) correlations in the 
nucleon momentum and removal energy spectra is adopted (see also 0,0). 

Theoretical predictions fall short of the data outside the interval 0.3 < a; < 0.55. It 
should be noted at this point that we have focused on the nucleon components only and 
that indeed Eq.(|l|) does not satisfy momentum conservation. This problem is generally 



^Note M\_^ = Ma-i + where E\_^ is the excitation energy of the outgoing A — \ nucleons 
system which is related to the nucleon's removal energy, E, hy E = Emm + Emm = Mjq + 

Ma-i - Ma. 

fA{z) properly includes the flux factor in its definition and it is normalized so as to satisfy baryon 
number conservation [^,11]. 



solved by introducing scattering from additional nuclear pions PJT^ which while restoring 
momentum conservation, also modifies the slope of the data at low x (at very low x nuclear 
shadowing eventually becomes dominant). Additional pions are not going to be of immediate 
relevance to our thesis which is to study the dependence of the moments of the non-singlet 
structure functions, for which the high x behavior is important. 

The next step is therefore to investigate whether the poor agreement with the data at 
high X is an indicator that lA is a poor approximation at > 2GeV^. In other words 
the effect of parton interactions can be sizeable and it can manifest itself differently inside a 
nucleus. We begin investigating this problem by defining the role of the transverse degrees 
of freedom and of the struck parton's virtuality (or off-shellness) in a nucleus. 

Notice that now both the struck parton and the parent nucleon are off-shell and the struck 
parton and spectator interactions can occur in more ways than in a free proton. With the 
intent of investigating this problem systematically and in a parton language suitable to 
nuclear targets, we start by comparing the on-shell and off-shell nucleon structure functions. 

The structure functions are extracted from the cross section (j(7*P —* k'X) {k' and X 
being the final quark and the spectator system). By using an infinite momentum frame 
parametrization of the particles four-momenta and Old-Fashioned-Perturbation-Theory 
(OFPT) one obtains 
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(2) 

where the struck parton's momentum components are = {xP + m^/2xP; px, xP), the 
spectators ones are Px = ((1 — x)P + m^/2(l — x)P; — px, (1 — x)P), a sum over the mx 
mass spectrum is implied in Eq.(^, 05 = S{x,p'^) = / x{l — x) +m\ / {1 — x) +171^ / x, (p{S) 
is the vertex function, and finally we are assuming scalar particles because our argument 
will be essentially kinematic (Eq.(^ can be easily generalized to spin 1/2 constituents). 
Notice also that the final result is Lorentz invariant and that the integration in Eq.(^) can 

^The mx distribution really determines the behavior of F2 only at x ~ anyway. 
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be directly transformed into an integral over S, thus transfering all the x dependence into 
the lower limit of integration: Smin = ^x/i^ — x) + rn? /x. 

By considering the same graph inside a nucleus, one obtains for the off-shell nucleon 
structure function 



l-v J 



y f pt 



MS) 



l-yj 2(27r)3 [y%S-M^^)\ " 
where now S = S{y, a^) = a^/[y{l—y)] +m\/ (1 — x) +rn? /y, with ax = Pt — ?/kT, y = x/ 
because of Fermi motion (the struck parton and the spectators have now perpendicular 
momenta: px and kx — Px? respectively, /ct being the nucleon transverse momentum). Two 
important observations emerge from the comparison of Eqs.(^ and (|^). First of all note 
the role played by the nucleon's kx'- while {z) ^ 1 and therefore rescaling effects in the 
longitudinal variables are predicted to be small, kx and px are similar in magnitude and we 
expect this to affect the derivation of the next to leading twist correlation functions in a 
nucleus |T^. Secondly one obtains {k^)/Mff modifications to the naive convolution formula 
(0), that can be calculated by relating F^^ and F^^^ by a shift of their longitudinal 
variables: F^^iy') = J{y,y')F2^^{y = x/z,kT), where y' y x {k^ + m\)/{yk'^ + m\), 
where J(?/', y) is a kinematical factor. This relation between -Pf ^ and F^^^ is obtained 
by transforming the integral over kx into an integral over S and by observing that the 
lower limit of integration, Smin = ^x/i^ ~ v) + Tn?/y + yVr/i^ ~ 2/)? is related to the 
free nucleon Smin-, calculated above by a shift of momentum and mass (this technique has 
been applied often in OFPT, see Ref. and references therein). It is valid within the 



hypothesis that the vertex function is not modified inside the nuclear medium, (0^ = 0), 
which is reasonable given the "locality" of the electron-nucleus interaction. The convolution 
integral that properly takes into account the transverse degrees of freedom (and the parton 
virtuality) is given by 

/•oo roo 

F^{x) = 2txMn dzz dkknA{k)F^^[y'{kl)l (4) 

Jx J k^in({E) A,z) 

where k = ^[MAr(l — z) — E]"^ + kx^, nA{k) is the nucleon momentum distribution []TT|, and 
both the Jacobian, 27r M^r, and the lower limit of integration, fc^m =| ^Af(l ~ ^) ~ {E)a \, 



come from the trasformation of variables (A;||,kT) — > {z,k). Eq.(^) satisfies, like Eq.(0), 
baryon number conservation but not momentum conservation because, as mentioned before, 
it only involves nucleons while it is well known that momentum can also be carried by nuclear 
pions. Eq.(|[) predicts a large EMC effect than Eq.(|I]) essentially because F^^{y') < {y), 
being y' > y. 

Therefore we conclude that in DIS scattering from a nucleus, the nucleon Fermi motion 
modifies the relationship between the struck parton virtuality and its transverse momentum, 
with respect to that of a free proton. It is well known that at high x the electron-nucleus 
cross section is sensitive to the high momentum components, i. e. to coherent effects arising 
from NN correlations that, as seen in Eqs.(0) and (|), can affect both the parton pt and x 
distributions. The physics at high x deals with the interplay between these coherent effects 
and the pure QCD interactions between the struck and spectator partons in the nucleus. This 
perhaps subtle point is the main result of the present paper which is intended to prepare the 
ground for our further investigation of higher twists. Our approach is substantially different 
from some commonly used ones that are based on a longitudinal convolution formula (Eq.(|I|)) 
while taking care of the off shell nucleon through ad hoc prescriptions. Such approaches are 
not only conceptually wrong, based on our previous discussion, but we show that they 
predict a pattern of x-scaling violations that is not observed in the data. In particular, we 
discuss the most popular model, obtained by extending to DIS the "De Forest's prescription" 
01 which was originally constructed for electron-nucleus Quasi-Elastic (QE) scattering and 
first applied to nuclear DIS in Ref. (see also for nuclear matter). The underlying 



philosophy is that lacking any knowledge on the form of the operators connecting off-shell 
nucleon states, one should reevaluate the kinematics according to: 



ko = MN-E- Ta-1 ^ko = + Ml, (5a) 

- {ko - ko) , (5b) 



where T4„i = A;^/2M4_i is the A — 1 system's recoil energy, ko is the off-shell nucleon 
energy whereas ko is DeForest's proposal for its on-shell substitute. This method essentially 
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reduces to putting the initial nucleon back on shell, at the expense of modifying the energy 
dependence of the four- momentum transfer, q^. Equivalently, the ambiguity due to the 
interaction involving an off-shell nucleon current is solved by "rescaling" the dependence 
of the nucleon structure functions. By implementing Eqs.(^, in the calculation of the 
inelastic cross section, one also obtains a breakdown of the convolution formula because of 



the modifications in the which is shared by the struck nucleon 

2. Since our final objective is to understand QCD in nuclei, in the last part of this 
paper we concentrate on the moments of the nuclear structure functions. It is a well known 
consequence of the lA with no off-shell effects that the nuclear moments factorize into 
the (Q^-dependent) nucleon moments and the (Q^-independent) moments of the nuclear 
light cone distribution. In what follows we show our theoretical expectations based on the 
previous discussion of off-shell effects along with our extraction of nuclear moments from 
current world experimental data. From our analysis we are able to determine the coefficient 
of the 0(1/(5^) correction and to compare it with the results of the quantitative analysis of 
Ref. IQ and which were obtained from data extending up to a lower value of x. 

The Cornwall-Norton moments of the nuclear structure functions read 

M^{Q^)= rdxx''~^F^{x,Q^). (6) 

JO 

If is written in terms of a longitudinal convolution formula the following factorization 
holds [ig 

' dzfAiz)z^-'\ 



dzfA{z)z''-' dxx''-^F^{x, = 
Jo Jo 



^M^M^iQ'), (7) 



which is such that the ratio M:^(Q )/M^(Q ) is independent of Q . A breakdown of the 
longitudinal convolution formula translates into a breakdown of such a factorization. We 
examined two different mechanisms deriving from off-shell effects, namely: (i) the presence 
of transverse degrees of freedom in the present approach, and (ii) the modification of the way 
Q"^ is shared by the struck nucleon according to Our results are presented in Fig.l for 



n = 3 — 5. All theoretical results were obtained by using the same realistic spectral function 
which includes NN correlations |T^. The experimental data are for ^^Fe; they have been 



obtained from an analysis of world data to be discussed later. One can see that while the 
convolution formula predicts no dependence at all, prescription (i) shows a very slight 
breaking of factorization which is evident from the small slope in and prescription (ii) 
seems to display an unphysical wiggle and a steeper slope, due to the way the is rescaled 
according to Eqs.(|5^^). The existing data are not accurate enough to detect any slope in 
the given range. 

In order to optimize the extraction of the dependence of we then investigated 
the following quantity: TZ^ = M^IQ"^) / M^{Ql) In the free nucleon case and for the 

non singlet (NS) part of the structure function, the LO perturbative QCD, prediction for 
this quantity is 



c^sm 



C (inQ'-lnAlcD) (8) 



asm 

where dn = •jn/'^Po, In are the NS anomalous dimensions, and Po = 11 — 23Nf. Ql is 
an arbitrary scale to be chosen conveniently. The motivation for normalizing the nuclear 
moments at a scale can be understood by looking at the behavior of the curves in Fig.l, 
where the different treatments of off-shell effects are shown to produce overall depletions 
of different size. In the ratio TZ^ such differences are leveled and one can concentrate on 
the Q^-dependence only. Having in mind Eq.(|^, in Fig. 2 we plotted TZ^ vs. InQ"^. For 
a free nucleon any deviation from a straight line can be considered as a signal of higher 
twist contributions, since it has been shown that NLO terms are negligible. For a bound 
nucleon one has the initial problem of distingushing the off-shell effects originating from 
the nucleon transverse degrees of freedom, from "pure QCD" higher twist effects. The 
data points shown are extracted from the results of experiments at CERN (NA2) |T8[ and 



SLAC (E139,NE18) p^,|20|. To combine these different sets, most of the data needed to be 



interpolated to common values of to generate measured curves. The NE18 data 
provided the high x values of F.f; however, the QE contribution was large and its removal by 



a calculation introduced additional systematic error. The extrapolation to values of x 
beyond the region of the data also introduced error. At low x, the error is small and becomes 
negligible for moments n > 2. At high x the extrapolation error is significant except where 
NE18 data is available. This error was estimated by looking at the range of values produced 
by different extrapolation methods. For all moments the extrapolation error is at least as 
large as the statistical uncertainty, and becomes dominant for large n. 

We compare data with our results for the convolution formula and our off-shell prescrip- 
tion. The normalization point is at = 12.5 GeV^. We present results for n > 3 because 
here the contribution of the singlet part of the structure function becomes negligible. The 
off-shell curve deviates very little from the free nucleon one Eq.(^ (they are indistinguish- 
able in the figure). Therefore we conclude that the off-shell effects considered here, i.e. the 
ones that derive from the modifications to the struck quark prp in a bound nucleon, have 
been properly subtracted through the use of the normalized ratios shown in Fig.2. 

We can now look for the presence of higher twist terms in the data. In order to address 
this problem quantitatively, we started from considering the simple parametrization 

M:^m = Mf^-^{Q')(^i + ^y (9) 

where M^^'~^^{Q'^) is given by Eq.(|D and is the mass parameter to be determined. We 
assume initially that the 0(1/(5^) correction is not A-dependent. Our results are presented 
in Fig.3 for a fixed range of values of Aqcd = 0.200 — 0.260 MeV. The values displayed show 
the sensitivity to the low region: results change dramatically by excluding the point at 
= 2 GeV^ (diamonds), and by including it (stars) where after n = 5 {i.e. where high x 
also dominates), the validity of Eq.(^) becomes doubtful. Our result is somewhat smaller 
than in the previous extraction, = (0.8 ± 0.2 GeV)^, 0, but it is consistent with the 
original estimate, = (0.4 GeV)^, and with Ref. [§] (see area between the dotted lines 
in Fig.3). 

In conclusion, we have shown that the complications of kinematical origin due to nuclear 
effects and nucleon off-shellness contribute a small, if not negligible, modification to the 



dependence of the moments of the nuclear structure function Therefore the consequences 
of QCD in nuclei can be examined in detail. In particular, we successfully extracted the 
coefficient of the 0(1/Q^) power corrections. As these effects are prominent in the region of 
low and large x, we anticipate the forthcoming Jefferson Laboratory experiments pl|j2^ 
to significantly reduce the experimental uncertainties in the extraction of the moments of 

We thank Ingo Sick for raising interesting questions. 
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FIGURES 

FIG. 1. Ratio of the moments for ^^Fe to the nucleon ones. The experimental data extracted 
in this paper are compared to: (a) the convolution formula; (b) the off-shell model described in 
this paper; (c) the De Forest prescription. 

FIG. 2. Pennington-Ross plots for the world data on ^2^. Solid curves: LO theoretical results 
using the convolution formula; The singlet contribution is shown for n = 3 (short dashed line) 
and it is negligible at higher n. The off-shell effects are shown for n=4 (short dashed line) as an 
expample since they are of the same magnitude, and therefore negligible, for the other values of n. 
Dotted curves: expectation for the "simple" higher twist contribution (Eq.(^)), calculated with 
from our experimental fit. 

FIG. 3. The mass parameter (squared), r^, which regulates the simple twist-4 expectation 
of Eq.(P). The diamonds are obtained from our fits to the 56Fe data, excluding the point at 

= 2GeV^; the stars include = 2GeV^; the two dotted lines show the range of the predictions 
by Ref. The value of originally obtained by de Rujula et al., was 0.16GeV^, which is still in 
qualitative agreement with our results. 
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